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Magnetization reversal by a femtosecond circularly polarized laser pulse has been 
recently demonstrated in rare-earth doped transition metals (RE-TM). The switch- 
ing mechanism has been attributed to an inverse Faraday effect and thermal effects. 
Based on the parameters provided in the experimental work, we show that this claim 
is unlikely to give rise to femtosecond reversal. Using a hybrid itinerant-localized pic- 
ture of the RE-TM system, we propose a new mechanism that requires the presence 
of the rare earth element to reduce the symmetry of the system as well as a strong 
enhancement of spin-orbit coupling between the d electrons and the f moments in 
the presence of the laser. Our model does not require the heating close to the Curie 
temperature of the sample. 

PACS numbers: 76.60.-hq,78.20.Ls,75.40.Gb, 76.60.Es, 52.38.-r 

Using perpendicular magnetized 6*^22-^674.66*03.4 thin films, it was shown that a 40 fs 
circularly polarized laser pulse with a wavelength of 800 nm is capable of switching the 
magnetization What is remarkable about this important and unexpected result is that 
it was all-optical and no external fields were needed. Instead the switching was critically 
dependent on the chirality and the intensity of the laser pulse. The minimum intensity of 
the laser needed to accomplish reversal was of the order of 10^^ W/cm^ which is still too low 
for the direct exchange of the angular momentum between the laser and the magnetization 
to be significant. The authors of Ref. [l| argued that a magnetic field generated by the 
inverse Faraday effect, HipE, in combination with the enhanced susceptibility at elevated 
temperatures is responsible for the switching. They estimated i^iPE to be of the order of 
10 T and the switching temperature Tsw to be close to the Curie point Tq. However an 

n 

equally valid reinterpretation of their data leads to Tsw ~ (^c ~ rroom)/2 ^]. The dynamics 



of the reversal was not investigated in but it is believed to occur on the femtosecond 
time scale. This belief is based on the fact that controlled switching is only possible if it is 
faster than relevant energy and phase relaxation times. 

In this paper, we first discuss the ideas put forward in Ref. l|] in a more quantitative 
way and show that their picture can not explain femtosecond magnetization reversal. After 
that, we suggest qualitatively a different richer picture of the dynamics in the system and 
point out the need to go beyond LLG to study magnetization dynamics on the femtosecond 
time scale. The necessity stems from the fact, that strong magnetic fields generated by 
the laser are alone not sufficient to drive femtosecond magnetization reversal. Instead it is 
the simultaneous optical enhancement of dissipation that enables ultrafast switching. This 
dissipation is made possible by the fast relaxation of the hot electrons which requires the 
full account of the non-adiabatic dynamics in the laser-spin-electrons system. This latter 
aspect of the problem is treated in detail in this paper since it is of tremendous importance 
to any non-equilibrium system and in particular fast spin dynamics in systems driven far 
from equilibrium by a laser pulse. Consequently, a self-consistent treatment of dissipation 
and fluctuations is required which is also included here. Our method of solution goes beyond 
the Born approximation and is non-perturbative which is needed to accurately treat spin 
dynamics in an interacting environment as discussed by DiVincenzo and Loss [3]. It is based 



on the two-point effective action 



4] and provides a systematic way to obtain solutions more 



accurate than those of the Bloch equations jsl. 

Before we discuss our model and method of solution, we first show that the ideas proposed 

n 

m Ref. which are based on LLG, can not form a basis for the understanding of the fast 
switching of the magnetization. The LLG equation is a phenomenological equation written 
for magnetization close to equilibrium and accounts for damping in the simplest possible 

way 

= -jMxWff + aMx . (1) 

dt dt ^ ^ 

The internal field and the external laser field, = Hi{cosuJit,smuJit,0), are part of H'^^^. 
In a frame rotating with the frequency ui = lO^^Hz, Hi becomes time independent and the 
lab frame time derivative ^ becomes ^ — t^/z x M. Hence in the rotating frame, the 
effective field is H^-'^-'^ -|- H^, where = —uJi/'-fz is the equivalent of the Barnett field which 
is of the order of 5 x 10^ Oe. This shows that the magnetization in the presence of the laser 
precesses at almost five orders of magnitude faster than in typical FMR experiments. The 
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LLG equation has been successful in reproducing most magnetization data with 'reasonable' 
damping constants a but it does a poor job in systematically accounting for damping in 
different correlated data sets. An example of such data is the measurement of damping in 
rare-earth doped transition metals (RE-TM) as a function of the orbital momentum of the 
impurities j^. 

The dynamics of the magnetization reversal as obtained from Eq.[T] with a = 0.2 and 
i^iFE = 10 T (experimental values of Ref. is shown in Fig. [1] for various laser fields. 
The laser intensity Hi/Hb = 0.001 corresponds to the experimental power and shows slow 
reversal within ^ 20 ps. The faster reversal curves for Hi/Hb = 0.01 and 0.1 is due to the 
Zeeman interaction with the laser but requires intensities well above the damage threshold. 
To achieve femtosecond switching at reasonable laser powers, the damping would have to be 
increased to non-physical values of a = 80.0 (cf. inset of figure [1]). This calls for a different 
way to include dissipation in fast spin dynamics. Heating GdFeCo to temperatures some 
tens of degrees below the Curie point also does not help achieve faster than picosecond 
reversal. The role of temperature is to increase the fluctuations in the initial values of Mz 
which will help speed up the reversal but not to the extent of being fast at the femtosecond 
scale. The Curve with initial magnetization Mz{0) = —0.5 shows the effect of increasing 
temperature on the switching speed. 

In the rest of this paper, we give a more careful treatment of the two aspects of time de- 
pendent fields and of damping in a simplified model accompanied by a qualitative discussion 
of the various factors involved in the reversal. 

The excitation of GdFeCo by the laser unquenches the orbital angular momentum of the 
hot electrons which are mainly d electrons provided by Fe. But the d electrons also carry 
magnetization and hence a strong interplay of orbital magnetization and spin magnetization 
in this system is expected. The Hamiltonian is 

H = Hd + Hph + Hi , (2) 
Hi = -eE{t) ■ r - H(t) ■ (/i^s + /i/S/) (3) 
+Js ■ S/ + Al ■ s + Xsol ■ Sf . 

The Hamiltonian Hd is that of the d (and s, p) electrons of both Fe and Gd, Hph is the 
Hamiltonian of the phonons and Hi is the interaction Hamiltonian which includes the dipole 
interaction, the Zeeman term with both d and f electrons, the spin-orbit coupling term of 
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FIG. 1: Dynamics of the perpendicular component in an inverse Faraday field Hjfe = 10 T. 
The curves are plotted for various laser fields Hi/Hb = 0.1, 0.01, and 0.001 with initial Mz{t = 
0) = —1. The curve strating at M^(0) = —0.5 emulates the effect of elevated temperatures on the 
reversal for Hi/Hb = 0.001. The LLG damping is a = 0.2. The inset is for Hi/Hb = 0.001 and 
the unphysical damping a = 80.0. 



the d electrons with coupling constant A, the antiferromagnetic coupling between the d and 
f electrons and finally the effective spin-orbit coupling between the orbit of the itinerant d 
electrons and the local f moment of Gd that gives rise to the perpendicular anisotropy in 
GdFeCo. This latter term has many common features to the well known Rashba coupling 
which gives rise to dissipationless spin currents in time-dependent or inhomogeneous fields 
[?]. It is the smallest of interactions in the absence of the laser but it lowers the symmetry of 
the effective field seen by the d electrons and stores the memory of the chirality of the light 
after the laser is turned off. The enormous complexity of the reversal process in GdFeCo is 
only partly reflected in this Hamiltonian and since no sufficient real-time data are available 
to warrant detailed calculations based on this model, we present here a calculation based on a 
further idealized Hamiltonian that still has the necessary elements needed for a femtosecond 
reversal process. Our reversal process is based on ideas similar to the (mechanical) Barnett 
effect jsl which says that an external torque applied to an iron rod affects its magnetization 
by acquiring a component along the axis of the torque. The torque in our case will be 
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provided by the electric field of the laser as it acts on the hot d electrons. However the 
magnetic field of the laser also provides a similar Barnett contribution through the time- 
dependent Zeeman interaction with the d and f moments. The phonons through the spin- 
orbit coupling and the dipole terms provide the dissipation required for the reversal process. 
Polarization in the presence of rotating fields and dissipation has lon g b een recognized in 
gravitational physics j^, accelerator physics 10 1 and in semiconductors 11 1. Therefore there 
is no surprise that our model will give rise to polarization induced by a circularly polarized 
laser. In order to give a clear discussion of the optical Barnett-like effect, we introduce the 
effective Hamiltonian, 7^^*^, for the spin degrees of freedom deduced from Eq.[2] and treat 
the orbital degrees of freedom and the effect of E; on them as part of the environment. 
TY*^// includes energy exchange between three different sub-systems (cf. Fig. [21 inset) and is 
{h=l) 

= ns + n, + Hs,i + nQ (4) 

The spin Hamiltonian, Tis = — 7H; ■ S — \ ASl, includes the interaction with and an axial 
anisotropy term. The S's represent the effective spin degrees of freedom in the GdFeCo 
system. The anisotropy, which is most relevant after the laser is turned off, is taken in the 
mean field approximation. The gyromagnetic ratio 7 is assumed positive. The Hamiltonian, 

2 2 

Tiq = ^ + ^q^, represents a single optical electronic mode with frequency ujq close to the 
aser frequency uji 12|]. The Hamiltonian Tiq is that of the macroscopic bath as defined in 
I3I . The interaction Hamiltonian Tiaqi = — AS-q — q-Q, includes an effective linear coupling 
of the spin to the optical mode through spin-orbit coupling, which is taken isotropic for sim- 
plicity, and a linear coupling of the optical mode to the phonons Q. The material dependent 
parameter A = Aq + \iEf / is an effective spin-orbit contribution with Aq ~ 0.001 eV, 
Ai ~ 0.1 eV and Ec is the critical field for switching which is taken from experiment. In the 
presence of the laser, an electron acquires additional orbital angular momentum e'^Ef /muf 
and hence the power dependence of A. As first pointed out in Ref. [l^, a strong spin-orbit 
coupling in NiO is needed to provide a fast relaxation channel for the spin of the electrons 
in the femtosecond regime. 

This effective Hamiltonian is the simplest model that describes well the femtosecond 
reversal and reproduces LLG close to equilibrium if c^o is much larger than the frequency 
of the spin [l^. In the present case, ujq is very close to that of the laser-driven spin and 
consequently one of the main assumptions underlying LLG, the separation of energy scales 
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FIG. 2: (a): Dynamics of the perpendicular component of the magnetization induced by 
circularly polarized light {Ei,Hi) with positive frequency ui = 10^^ Hz and various intensities 
Hi/Hb = 0.00001, 0.0005, 0.001, and 0.002. The initial magnetization is M^O) = -1. (b): The 
dependence of switching on the chirality and initial conditions = ±1 for Hi/Hb = 0.001. All 
results are for ujQ = 0.8 uou r = 0.2, and A = 10^ Oe. Hi = IO'^Hb corresponds to a power of 
^ 5 • lO^^W/cm^. 

between the spin and the environment, is broken [l^. Hence we expect non-adiabatic 
behavior to be important in this system. We will show that a careful treatment of the 
dynamics give non-exponential switching different than the one expected from the Bloch 
equations [5| or the LLG equation. This is the physics we would like to address in our 
analysis of the switching since similar conditions are also present in other important problems 
such as those in quantum computation where it may be necessary to go beyond the Born 
approximation to study coherence j^. 

The coupled spin-laser-bath system is better studied in a frame (xi, X2, X3 = z) rotating 
around the 2;— axis with frequency uJi. In this frame, the Larmor torque is time-independent 
and the spin Hamiltonian becomes 



H' 



■^Mi - Xi:,MA (t) -^{Bo + Hb) M3 



(5) 



where Bq = A (Sz). The equations of motion for (M), the average of the spin operator, are 

(Ml) = -7(5o + Hb){M2) + 2X{{M2h) - (M362)) (6) 
(M2) = 7(^0 + Hb){M^) - ^Hi{M,) + 2X{{M,h) - (M^b,)) 
(Ma) = ^Hi{M2) - 2X ((M261) - (M162)) , 
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where 6j's are the harmonic mode variables in the rotating frame. The equations of motion 
for q are 

(^+a;o^)(g,) = (Q,) + A(cr,), (7) 

The bath Q is assumed Ohmic and is the source of dissipation in the mode q which we 
take to be r = 0.2. In the adiabatic hmit and in the absence of anisotropy, the system spin 
plus bath, Hs + Hq, has been treated earlier, see e.g. Ref. [l^]. Here we study the more 
complex and experimentally relevant case of non-adiabatic magnetization reversal where the 
inclusion of memory effects is necessary. 

In the rotating frame, the instantaneous torque is modified by the Barnett field. For 
visible light with ui ^ lO^^Hz, the optical Barnett field is as large as ~ 10^ Oe. This is 
larger than most exchange fields and hence a single spin picture should be adequate even for 
the treatment of the interaction of a laser with a ferromagnet. Despite its tremendous mag- 
nitude, we find that this effective field does not induce femtosecond magnetization reversal 
unless three key requirements are met. First, the Barnett field has to be much larger than 
the laser field H^. Secondly, the magnetization has to be coupled through the electric field 
E; to at least one optical mode q of energy c^o ~ ^Hb- Third, the damping of this mode due 
to its interaction with the macroscopic bath Q has to enable efficient energy transfer from 
the spins to the mode q. All three requirements can be met by various combinations of the 
parameters, Hb, Hi,Ei, uq and A, but the range of each individual parameter depends on 
the particular choice of the others. 

To solve equations [6] and U\ we need to calculate the average of the product of two 
operators {ai{t)qj{t')) . For this we need the density matrix of the whole system or we may 



use the more useful functional formalism instead 



for a detailed discussion of this method 



and its application to sd exchange in metals see 18|] and references therein). The generating 
functional is 

Z[Ji,J2] = (^j Dr]DpD(iexp (^-i j dtiHsiv^) (8) 
(p, q) + Hsqi - Ji ■ q - J2 ■ cr)))^ 

where the spin variables are for simplicity assumed 1/2 and are expressed in terms of Grass- 
mann variables, a = —^rjxr] [l9|. This allows the use of Wick's theorem in the path integral 
expansion and Ji are two virtual external sources. The problem then becomes very similar 
to that of a Fermi gas which was discussed with the aid of two point functions {4]. The 
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average values are calculated in the same way as in Ref. ^ 

Before switching on the laser, the phonons and the magnetization are in equilibrium and are 
weakly coupled since the electrons are in their ground state with a quenched orbital angular 
momentum. After that, the spin is driven out of equilibrium by the laser and q is treated 
as a perturbation which responds linearly to any changes in the spin. As in the sd exchange 
problem [l^, the equations of motion are non-local in time and integrating out the optical 
mode q gives rise to dissipation and fluctuations even at zero temperature. 

Using a Runge-Kutta scheme, Eqs.[6] and [7] are solved self-consistently. In line with 
the discussion above, we find that switching occurs only for a limited range of mutually 
dependent parameters. The results in fig.[2] focus on the power dependence of the magnetic 
response for uji = 10^^ Hz, uq = O.Slji, T = 0.2, and A = 10^ Oe. While the response to laser 
fields Hi/Hb < 10~^ is negligible on a femtosecond time scale, it becomes significant for 
Hi/ Hb = 5 • 10^'* and ultrafast reversal is found at higher fields of 1Q~^Hb and 2 • lO^^if^. 
At the highest fields, shows an almost instantaneous reversal accompanied by strong 
oscillations. The period of these oscillations is governed mostly by Hb while their decay 
depends on F. The reversal slows down with time and the approach to equilibrium depends 
on the coupling A of the spins to the mode q which is proportional to the power. The curves 
in the inset show the strong dependence of the reversal on the chirality of the laser. It is also 
important to realize that the simple exponential behavior expected from the Bloch equation 
or the LLG equation for the decay of the z-component is no longer true in our model. Fig. 
[2] shows that exponential behavior is only recovered after the laser has been on for several 
periods. The initial response of the spin system is therefore highly non-adiabatic due to 
the simultaneous action of the laser and of the dissipation and fluctuations, which are also 
enhanced by the laser, due to the coupling of the hot electrons to the magnetization. 

In summary, we have shown that circularly polarized light can induce femtosecond magne- 
tization reversal, when optical electron-phonon modes with frequencies comparable to those 
of the light are made available by the electric field of the laser through spin-orbit coupling. 
We predict a fast non-exponential reversal that cannot be recovered from the modified Bloch 
equations or the LLG equations and should be observable in future time-resolved experi- 
ments. The switching does not require high temperatures and is insentitive to the anisotropy 
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which makes the proposed mechanism very attractive to high density magnetic recording. 
We thank D. Stanciu, A. Kirilyuk and Th. Rasing for discussing their results with us. 
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